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Abstract 

We study boundary value problems posed in a semistrip for the elliptic sine-Gordon 
equation, which is the paradigm of an elliptic integrable PDE in two variables. We use 
the method introduced by one of the authors, which provides a substantial generaliza- 
tion of the inverse scattering transform and can be used for the analysis of boundary 
^ I as opposed to initial-value problems. We first express the solution in terms of a 2 x 2 

. matrix Riemann-Hilbert problem formulated in terms of both the Dirichlet and the 

Neumann boundary values on the boundary of a semistrip. We then concentrate on 
the case that the prescribed boundary conditions are zero along the unbounded sides of 
the semistrip and constant along the bounded side; in this particular case we show that 
CN ■ the "jump matrices" of the above Riemann-Hilbert problem can be expressed explicitly 

in terms of the width of the semistrip and the constant value of the solution along the 
bounded side. This Riemann-Hilbert problem has a unique solution. 



1 Introduction 

A method for solving initial-boundary value problems for linear and integrable nonlinear 
PDEs was introduced in 13J and developed by several authors T^. This method has 
already been used for: 

(a) linear and integrable nonlinear evolution PDEs formulated on the half line and on a 
finite interval [1 El [III [11 [11 [Ti [111 [13 [M 113 ISSl ISZl [Ml SS] ; 

(b) linear and integrable nonlinear hyperbolic PDEs [26l [41] ; 

(c) linear elliptic PDEs [l[Hl[a[Iil[Ml[5Hl[3nill3lllll5]. 

The aim of this paper is to implement this method in the case of the prototypical integrable 
nonlinear elliptic PDE, namely the celebrated sine-Gordon equation ( simple boundary value 
problems for this equation, using the method of 13 , have been analyzed in [TO I42j ). 
We will consider the sine-Gordon equation in the form 

Qxx + Qyv ^ s'mq, q^q{x,y), (1.1) 
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Figure 1: The semistrip S 

and we will analyze boundary value problems posed in the semi-infinite strip 

iS = {0 < a: < oo, <y < L}, 

where L is a positive finite constant. The sides {y = L, < x < oo}, {a; = 0, < y < L} 
and {y = 0, < x < oo} will be referred to as side (1), (2) and (3) respectively, see figure 1. 
Suppose that (|l.ip is supplemented with appropriate boundary conditions on the boundary 
of the semistrip S, so that there exists a unique solution q(x, y). It will be shown in section 2 
that this solution can be expressed in terms of the solution of a 2 x 2 matrix Riemann-Hilbert 
(RH) problem with jumps on the union of the real and imaginary axis of the A complex 
plane. The "jump matrices" are expressed in terms of certain functions, called spectral 
functions, which will be denoted by {aj(A), 6j (A)}, j = 1,2,3. These functions can be 
uniquely characterized via the solution of certain linear Volterra integral equations, in terms 
of the Dirichlet and Neumann boundary values. Namely, {ai,bi}, {02,62} and {03,63} are 
uniquely determined in terms of {q{x, L), qy{x, L)}, {q{0,y), qx{0,y)} and {q(x, 0), (7j,(x, 0)} 
respectively. However, for a well posed problem only a subset of these boundary values are 
prescribed as boundary conditions. Thus, in order to compute the spectral functions in terms 
of the given boundary conditions, one must first determine the unknown boundary values, 
i.e. one must characterize the Dirichlet to Neumann map. The solution of this problem, 
which makes crucial use of the so-called global relation, yields in general a nonlinear map. 
In the case of integrable nonlinear evolution PDEs, it has been shown in [TBI [HI EH [35] 
that there exists a particular class of boundary conditions, called linearizable, for which it 
is possible to avoid the above nonlinear map. The main result of the present paper is the 
analysis of linearizable boundary conditions for the sine-Gordon equation on the semi- infinite 
strip. In particular, the following boundary conditions will be investigated in detail: 

q{x,L) = q{x,0) = 0, < x < 00; q{0,y)=d, 0<y<L, (1.2) 

where d is a finite constant. It will be shown in section 5 that in this particular case, the 
"jump matrices" of the associated RH problem can be constructed explicitly in terms of the 
given constant d and of the width L of the semistrip. This result, as well as the analogous 
result valid for the elliptic version of the Ernst equation [25] , imply that the new method of 
[13] provides a powerful tool for analyzing effectively a large class of interesting boundary 
conditions. 
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2 Spectral analysis under the assumption of existence 

In what follows we assume that (ll.ip is supplemented with appropriate boundary conditions 
on the boundary of the semistrip iS so that the existence of a unique solution q{x,y) can be 
assumed. Furthermore, we assume the following: 

q{x,L), qy{x,L), q{x,0), qy {x, 0) e h\R+), 

xq{x,L), xqy{x,L), xq{x,0), xqy{x,0) £ L^{S.'^), (2-1) 
q{0,y), g,(0,y), yqiO,y), yfe(0,j/) G L\[0,L]). 

The sine-Gordon equation is the compatibility condition of the following Lax pair for the 
2x2 matrix-valued function '^{x, y, A), A G C: 

n(X) 

*x + -^[^3,*] =Q(a;,y,A)vI/, (2.2) 
^y + ^[a3,^]^iQix,y,-X)^, (2.3) 

where 

1 ^, 1\ 1 1 



^^(A) = -^A--J, .(A) = -^A + -J, (2.4) 

^/ i(l-cosg) + 

Q(a;,?/,A) = - , q = q{x,y). (2.5) 
\qx-iqy-'-^ -i(l-cosg) / 

Equations (|2.2p and (|2.3|) can be written as the single equation 

d (^e(^(^)^+'"(^)2')^) ^Sj{x, y, A) = Wix, y, A), (2.6) 

where the differential form W is given by 

W{x, y, A) = e(^^(^)^+"(^)^)^' {Q{x, y, X)'f{x, y, X)dx + iQ{x, y, -X)'f{x, y, X)dy) , (2.7) 
and CT3 acts on a 2 x 2 matrix A by 

[a3,A]. 



Remark 2.1 Note that 



n{X) = -n{X) = n{-), a;(A) = w(A) = t^(-). 
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Figure 2: The functions ^'^ 



2.1 Bounded and analytic eigenfunctions 

We define three solutions '^j{x, y, X) j — 1,2, 3, of (|2.6p by 

■^,ix,y,X) = I+ /"^"'''\-("(^)"+"(^)^)^VF(C,ry,A), (2.8) 

where 

= (oo,y), (x2, ya) = (0, L), (xg, ya) = (0, 0). (2.9) 

Since the differential form W is exact, the integral on the right hand side of (|2.8p is inde- 
pendent of the path of integration. We choose the particular contours shown in figure 2. 
This choice implies the following inequalities on the contours: 

{xi,yi) ^ {x,y) : ^ - a; > 0. 

{x2,y2) ^ {x,y) : X <{), ri-y>Q. 

{X3,y3) ^ {x,y) ■■ ^ ~ X <0, ri -y <0. 

The first inequality above implies that the exponential appearing in the second (first) column 
of the right hand side of the equation defining is bounded and analytic for Im{X) < 
(Jm(A) > 0). Similar considerations are vahd for ^'2 and ^3. Hence we denote the matrices 
as follows: 

vE'i = (*(^'\*f^ *2 = (*^'\*f), *3 = (*f 

where the superscript (12) denotes the union of the first and second quadrants of the A 
complex plane, and similarly for the other superscripts. The function ' is analytic for 
Im{X) > and it has essential singularities dA, X — 00 and A = 0; furthermore, 

( J )+o(^), A ^00, /m(A)>0. (2.10) 
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Similar considerations are valid for the column vectors ^'g'^-' and 5*3^^ The function 

^2 is an analytic function in the entire complex plane, except at A = cxd and A = 0, where 
it has essential singularities. In addition, 

*^'^-( A^oo, ^<arg{X)<2n, 

M-f = J ^+oQ), A^^, ^<arg{X)<7r. (2.11) 

2.2 Spectral functions 

Any two solutions of (|2.6p are related by an equation of the form 

y, A) = ^{x, y, X)e-^"^^^^+^^^^y^^ C{X). (2.12) 

We introduce the notations 

^i(A) = *i(0,L,A), ^2(A) =*2(0,0,A), ^3(A) = «'i(0,0,A). (2.13) 
Then equation (|2.12p implies the following equations: 

*i(a;,y,A) = «'2(x,j/,A)e-(^^(^'^+"(^'^)^e'^^^S'i(A), Ae (R+,R"), (2.14) 
^'2(x,y,A) = *3(x,2/,A)e-(^^(^)"+'^(^)^)^'52(A), Ae(^M-,^M+), (2.15) 
M'i(a;,y,A) = vl/3(x,y,A)e-("(^)"+'^(^>y)*53(A), Ae(R-,M+). (2.16) 

The notation A S (M^ , ) means that the equation for first column vector of (|2.14p is valid 
for A € R"*", while the equation for the second vector is valid for R~, and similarly for (|2.15p . 

Equations (|2.13p - (|2.16p suggest the following definitions: 

poo 

^i(A) = <i>i(0,A), $i(x,A) = /-/ c^^(^)(«--)^g(e,L,A)$i(e,A)de, 

J X 

Ae(C+,C"), 0<x<oo, (2.17) 

rL 

S2iX) = <i>2(0, A), <i>2(y, X) = I-i e"(^)(''-^)^g(0, 77, -X)^2{v, >^)dv, 

AgC, 0<y<L, (2.18) 

^3(A) = ci>3(0,A), ^3{x,X)=I- e"(^)(«-^)^Q(e,0,A)a>3(C,A)dC, 

J X 

Ae(C+,C"), Q<x<oo. (2.19) 
The matrix Q satisfies the symmetry properties 

Q(A)22 = Q(-A)ii, Q(A)i2 = Q(-A)2i. 
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Hence the matrices $i can be represented in the form 



f A,{x,X) Bi{x,-X)\ ( A2{y,X) B2{y,~X) \ {A3{x,X) B3{x,~X) \ 

\Biix,X) Ai{x,-X) J' ^2 \B2iy,X) A2iy,-X) J ' \B3ix,X) A3{x,-X)J 



The spectral functions {ai(A), 61(A)}, {a2(A), 62(A)} and {a3(A), 63(A)} are defined in 
terms of { g(x, £ ) , qy{x ,L)}, {q{0,y), qx{0,y)} and {q{x,0), qy{x,0)} respectively, through 
equations ^JT^-^J^. 

These functions have the following properties: 

• ai(A), 61(A) are analytic and bounded in C^. 
ai(A)ai(-A) -6i(A)6i(-A) = 1. 

ai(A) = 1 + O (i), 61(A) = O (i) as A ^ 00, Im{X) > 0. 

• a2(A), 62(A) are analytic functions of A for all A G C, except for essential singularities 
at A = 00 and A = 0. 

a2(A)a2(-A) -62(A)62(-A) = 1. 

a2(A) = 1 + O (i), 62(A) = O (i) as A ^ 00, ^ < arg{X) < 2tt. 

• 03 (A), 63(A) are analytic and bounded in C^. 
a3(A)a3(-A)-63(A)63(-A) = l. 

a3(A) = 1 + O (i), 63(A) = O (i) as A ^ 00, Im{X) > 0. 

These properties follow from the analogous properties of the matrix- valued functions 

j — 1,2,3, from the condition of unit determinant, and from the large A asymptotics of 

these functions. 

2.3 The global relation 

Evaluating equations (|2.15p and (|2.16|) at x = 0, y = L, wc find 



and therefore 




i = 1,2,3. 



J=*3(0,L,A)e~^^"^52(A) 



and 



Si{X) = *3(0,L,A)e-^^^^ 53(A). 



Eliminating ^'3(0, L, A) we obtain 



SiiX)^S2iX)-^S3{X). 



(2.20) 



The first column vector of this equation yields the following global relations: 



ai(A)=a2(-A)a3(A)-62(-A)63(A), 
6i(A)e--(^)^ - a2(A)63(A) - a3(A)62(A), 



A e C+, 
A e C+. 



(2.21) 
(2.22) 
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2.4 The Riemann-Hilbert problem 

Equations (|2.14p - (|2.16p . relating the various analytic eigenfunctions, can be rewritten in a 
form that determines the jump conditions of a 2 x 2 RH problem, with unitary jump ma- 
trices on the real and imaginary axes. This involves tedious but straightforward algebraic 
manipulations. The final form is 

'i'^{x,y,X) = '^+{x,y,X)J{x,y,X),X eRliiR, (2.23) 

where the matrices and J are defined as follows: 

,(12) 1 ,T,(2) \ ... 



2 

J{x,y,X) ^ J°'{x,y,X), if arg{X) a, a = 0, |, tt, (2.24) 



where, using the global relation, we find 



( 



a2W baj-X) ^-g(x,y.X) 

ai(-A)a3(A) a3(A) 

j°= : 

\ ai(-A) 



h2(-A) ^-e(x.t,.A) 



fc2(A) ^e(x.v.\) I 
ai( — A)a3( — A) 



and 

where 



9{x,y,X) ^n{X)x + uj{X)y. (2.26) 

All the matrices J" have unit determinant: for J^/^ and J^'^/^ this is immediate, whereas 
for J*^ we find 

jg^(^O) ^ a2(A)+e-"(^)^6i(A)63(-A) ^ ai(-A)a3(A) ^ 
ai(-A)a3(A) ai(-A)a3(A) 

where we have used the equation 

a2(A) =ai(-A)a3(A)-63(-A)6i(A)e-"(^)^, A £ M. (2.27) 
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Figure 3: Bounded eigenfunctions and the Ricmann-Hilbert problem 



Equation (|2.27p is a consequence of equations (|2.2ip and (|2.22p (see also equation (|4.19p 
below) . 

The function ^{x,y,X) solution of this RH problem is a sectionally meromorphic function 
of A. The possible poles of this function are generated by the zeros of the function ai(A) in 
the region {arg{X) € [^,7r]}, by the zeros of a3(A) in the region {arg{X) G [0, ^]}, and by 
the corresponding zeros of ai(— A), a3(— A). 

We assume 

• The possible zeros of ai in the region {arg{X) G [f,7r]} are simple; these zeros are 
denoted Aj, j = 1, .., A^i 

(2.28) 

• The possible zeros of 03 in the region {arg{X) G [0, 1]} are simple; these zeros are 
denoted Q, j — 1, ..,N3 

The residues of the function ^! at the corresponding poles can be computed using equations 
((m)) - ((TTB)) . Indeed, equation ((TTB)) yields 

hence 

a3 - d3(0) "a3(0)63(0) ' ^ ^ ^ 

where 03 (A) denotes the derivative of 03 with respect to A. 
Similarly, using (|2.14|) . 

ai " ai(A,) ~ ai(A,)6i(A,)e-(^.)i ' ^ ^ ^ 

Using the notation for the first column, [5']2 for the second column for the solution 
^' of the RH problem (|2.23p . at equations (|2.29p and (|2.30p imply the following residue 



8 



conditions: 



(2.31) 



and similar ones in C by letting A ^ — A. 

The inverse problem 

Rewriting the jump condition, wc obtain 

- = *+ - J- = «'+(/- J) ^ - = (2.32) 
where J = / — J. The asymptotic conditions (|2.10p - (|2.1ip ) imply 

^i,{x,y,X) = I+^-^ + o(^), \X\^oo. (2.33) 



A VA^, 

Equations (|2.32p and (|2.33p define a Riemann-Hilbert problem. 
The solution of this RH problem is given by 



^(x,y,A)^/ + ^/^ ^+(--'^;yf Aer, (2.34) 



/r 
where 

r = MU; 

Equations ((OS)) and fOil imply 



= -T^ / ^+ix,y,X)Jix,y,X)dX. (2.35) 



r 



27ri 

Using (p:^ in the first ODE in the Lax pair (|^ . we find 



- ^[(73, - z^^-^ai ^ - zg, = 2(**)2i = 2 lim (A*2i), (2.36) 

4 4 A— ►oo 

(ci, CTs denote the usual Pauli matrices). 

In order to obtain an expression in terms of q rather than its derivatives, we consider the 
coefficient of the term A~^. The (1,1) element of this coefficient yields 

cosg(a;, y) = l + MK)n - 2(v|/*)2^. (2.37) 

3 Spectral theory assuming the vahdity of the global 
relation 

3.1 The spectral functions 

The above analysis motivates the following definitions for the spectral functions. 
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The spectral functions at the y = and y = L boundaries 

Definition 3.1 Given the functions q{x,L), qy{x,L) satisfying conditions i2.1]) . define the 
map 

Si : {q{x,L),qy{x,L)} {ai(A), 5i(A)} 

by 

' ai(A) 



61(A) 



[$i(0,L)]i, Ae 



where [$i(x, L)]i denotes the first column vector of the unique solution $i(x, L) of the 
Volterra linear integral equation 

roc 

$(x,A)=/-/ e"(^)(«--)^Q(^,L,A)ci>(e,A)dC, (3.1) 



and Q{x, L, X) is given in terms ofq{x,L) and qy(x, L) by equation I12.5\) . 
Proposition 3.1 The spectral functions ai(A), 61(A) have the following properties. 

(i) ai(A), fei(A) are continuous and bounded for I m{\) > 0, and analytic for I m{X > 0. 

(li) ai(A) = 1 + O (i), fei(A) = O (j) as A ^ 00, /m(A) > 0. 

(iii) ai(A) = cos + 0(A), 61(A) = i sin + O (A) as A -> 0, Im{X) > 0. 

fi?;; ai(A)ai(-A) - 6i(A)6i(-A) = 1, /m(A) > 0. 

(v) The map Qi : {oi, 61} ^ {q{x, L) qy{x, L)}, inverse to Si, is given by 
cosq{x, L) = 1 + 4j lim (AMa;)ii + 2 lim (AM)2i, 

A^oo A— >oo 

qy[x,L) = -iq^{x,L) + 2 lim (AAf)2i, 

A — >oo 

where M is the solution of the following Riemann-Hilbert problem: 

* The function 

is a sectionally meromorphic function of X ^ <C. 

* M ^I + 0{{) as X-^ 00, and 

M+{x,X) = M^{x,X)Ji{x,X}, AeR, 

where 

^i(^,A)=K^ , AGM. (3.2) 

V ai(-A)^ ai(A)ai(-A) / 

* The function ai(A) may have Ni simple poles Xj in C+. 
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* Let [M]i denote the i-th column vector of M , 1 — 1,2. The possible poles of 
occur at Xj, and the possible poles of M- occur at — Aj in C~, and the associated 
residues are given by 

Resx^[M{x,\)]2 = . \ [M(x,\j)]i, 
ai(Aj)bi(Aj) 

i?es_A,[M(a:,A)]i = ■ ( f ^ X J M(x, -A,-)]2. (3.3) 
The spectral functions {03,63} satisfy an analogous result: 

Definition 3.2 Given the functions q{x,0), qy{x,0), satisfying conditions i2.1\) . define the 
map 

§3 : {g(a;,0),g,(a;,0)}->{a3(A), 63(A)} 

by 

( 53(A) ) = [*3(0,0)]„ AeC+, 

where [$3(0;, 0)]i denotes the first column vector of the unique solution $3(0;, 0) of the 
Volterra linear integral equation 

/•oo _ 

$(x,A)=/-/ e"(^)(«-^)*Q(e,0,A)$(e,A)de, (3.4) 

and Q{x,0, X) is given in terms ofq{x,0) and qy{x,0) by equation i2.5\) . 

Proposition 3.2 The spectral functions a3(A), 63(A) have the properties (i)-(v) of proposi- 
tion iS. 1]) . provided ai is replaced by as, 61 is replaced by 63, Si is replaced by S3 and L in 
replaced by in all expressions. 

The spectral functions at the x = Q boundary 

Definition 3.3 Given the functions q{0,y), qx{0,y), satisfying conditions i2.1]) . define the 
map 

S2:{g(0,y),g,(0,y)}^{a2(A), 62(A)} 

by 

( 62(A) )=t*^^°'"^^^' 

where [$2(0,y)]i denotes the first column vector of the unique solution $2(0, y) of the 
Volterra linear integral equation 

Hy, X)=I-if e-(^)(''-^)*Q(0, y, -A)<i>(7y, X)d7^, (3.5) 

and Q{0,y,X) is given in terms of q{0,y) and qx{0,y) by equation i2.5\) . 
Proposition 3.3 The spectral functions a2(A), 62(A) have the following properties. 
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(i) a2(A), 52(A) are entire functions of X, except for essential singularities at X = and 
X ~ 00, bounded for Re{X) > 0. 

(li) a2(A) = 1 + O (i), b2{X) = O (j) as X^ 00, Re{X) > 0. 

(iii) a^iX) = cos + O (A), 62(A) = i sin + O (A) as A ^ 0, /m(A) > 0. 

(iv) a2(A)a2(-A) - 62(A)52(-A) = 1, /m(A) > 0. 

(^■yj The map Q2 ■ {02, 62} ^ {9(0, y) 9i/(0, y)}, inverse to §2, is given by 
cosq{0,y) = I + Ai lim iXMy)n +2 lim (AM)2i, 

A— >00 A— >CXD 

g:,(0,y) =ig,j(0,z/) + 2 lim (AAf)2i, 

A — »oo 

where Ai is the solution of the following Riemann-Hilhert problem: 

* The function 

M(vX)^{ ^^^^0 
^vnv,^) <y M^{y^X) ReX<0 

is a sectionally meromorphic function 0/ A G C. 

* M = 1 + (i) as A ^ 00, and 

M+(y,A) = Af_(2/,A)J2(y,A), X e iR, 

where 

( 1 fc2(-A) ^-^a):r \ 

^2(y,A)=. ^1 , AGzM. 

\ a2(-A)^ a2(A)a2(-A) / 

* M satisfies appropriate residue conditions at the zeros of a2{X). 
Proof of propositions (f3TT|) - (f3T3l) 

The proof of properties (i)-(iv) follows from the discussion in Section 2.2. In particular, 
property (iii) follows from the asymptotic behaviour at A ^ 0, which can be derived by 
analysing equations (|2.2p - (|2.3p (see [ID]), and is given by 

vI/ = vI/o + 0(A), |AHO, *o(:r,y)= ( V (3-6) 

\ I sm cos J 

To prove (v), we note that the function (t>i{x, A) is the unique solution of the ODE 

(kx + ^^^'Slcj) = Q(x, L, A)(/)(x, A), 
lim 4'{x, A) = /. 

A—* 00 

Furthermore, 4>-i{x, A) is the solution of the same ODE problem, with Q{x, L, A) replaced by 
Q(a;,0,A). 
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Similarly, <j>2{x, A) is the unique solution of the ODE 



(j)y + CT3(/> = iQ{0, y, -A)(/i(x, A), 
</)(0,L)=/. 

The spectral analysis of the above ODEs yields the desired result. 

Regarding the rigorous derivation of the above results, we note the following: If {q{x, L),qy{x, L)}, 
{q{x,0),qy{x,0)} and {q{y,0),qx{y:0)} are inL^, then the Volterra integral equations p.ip . 
(|3.4p and (|3.5p respectively, have a unique solution, and hence the spectral functions {aj, bj}, 
j = 1,..,3, are well defined. Moreover, under the assumption (|2.ip the spectral functions 
belong to H^(R) , hence the Riemann-Hilbert problems that determine the inverse maps 
can be characterized through the solutions of a Frcdholm integral equation, see [TUl I49j . 
QED 



3.2 The Riemann-Hilbert problem 

Theorem 3.1 Suppose that a subset of the boundary values {q{x, L), qy{x, L)}, {q{x, 0), qy{x, 0)}, 
< a; < oo, and {q{y,0),qx{y,0)}, < y < L, satisfying are prescribed as bound- 

ary conditions. Suppose that these prescribed boundary conditions are such that the global 
relations \2.21\) and \2. 22\) can be used to characterize the remaining boundary values. 
Define the spectral functions {aj,bj}, j = 1, ..,3, by definitions i3.1\) - l[KM) . Assume that the 
possible zeros o/ai(A) and {Cj}f2i o/a3(A) are as in assumvtion \2.2^ 

Define M{x, y, A) as the solution of the following 2x2 matrix Riemann-Hilbert problem: 

* The function M (x, A) is a sectionally meromorphic function of A away from M U iM. 

* The possible poles of the second column of M occur at \ ~ Q , j ^ 1, N2, in the first 
quadrant and at X = Xj , j ~ 1, Ni, in the second quadrant of the complex X plane. 

The possible poles of the first column of M occur at X = —Xj (j = l,...,Ni) and 
X^-Q (j^l,...,N2). 

The associated residue conditions satisfy the relations 12.31]) . 

* M = I + 0{{) as X^ 00, and 

M+{x,y,X) = M^{x,y,X)J{x,y,X), A e R U 

where M = for X in the first or third quadrant, and M = Af_ for X in the second 
or fourth quadrant of the complex X plane, and J is defined in terms of {aj,bj} by 
equations \2.25\) . 

Then M exists and is unique, provided that the H"'^ norm of of the spectral functions is 
sufficiently small. 

Define q{x,y) is terms ofM(x,y,X) by 

qx-iqy ^ 2 lim (AA/)2i, (3.7) 

A — ^OQ 

cosq{x,y) = 1 + 4i( lim (AA-f:r)ii) - 2( lim (AAf)2i)^ (3.8) 

A — >oo A — >oQ 
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Then q(x,y) solves Hl.lfl . Furthermore, q{x,y) evaluated at the boundary, yields the func- 
tions used for the computation of the spectral functions. 

Proof: Under the assumptions (|2.ip . the spectral functions are in H^. 
In the case when ai(A) and 03 (A) have no zeros, the Riemann-Hilbert problem is regular 
and it is equivalent to a Fredholm integral equation. However, we have not been able to 
establish a vanishing lemma, hence we require a small norm assumption for solvability. 
If ai(A) and 03 (A) have zeros, the singular RH problem can be mapped to a regular one 
coupled with a system of algebraic equations [22]. Moreover, it follows from standard ar- 
guments, using the dressing method [47l [48] , that if M solves the above RH problem and 
q{x,y) is defined by (|3.7p - (|3.8p . then q{x,y) solves equation The proof that q evalu- 

ated at the boundary yields the functions used for the computation of the spectral functions 
follows arguments similar to the ones used in [24j . 
QED 



d^, 



4 Linearizable boundary conditions 

We now concentrate on the particular boundary conditions (|1.2p . 
In this case, equations (|2.17p - (|2.19p simplify as follows: 

Aiix,X) \ ^ f r { g,(e,L)Ai(e,A) 

< X < 00, Im{\) > 0, 
A^{y,\) \ ^ /^iVi/V -^^^A^{^,X) + [q,{Q,y)~^^^]B^{^,\) 
B2{y,\) ) ) ^Jy \ e-(^)(^-'') [fe(0,y) + i2iaii]A2(^,A)+ (i^i?2(r7,A) 

Q<y <L, A e C, 

^3(x,A) \ _ (^\_\r( 'Zy(e,0M3(e,A) 
Bz{x,\) ) \0 J aJ, [ e"W(-«)g,(e,0)i?3(e,A) 

< x < 00, /m(A) > 0. 

In equations (|4.ip and (14. 3p . the only dependence on A is through il(A). Thus, since f2( — j) — 
J7(A), it follows that the vector fmictions {Ai,Bi) and (A^jB^) satisfy the same symmetry 
properties. Hence, 

a,(-i) =a,(A), 6^.(-l) = 6^-(A), i = l,3, /m(A) > 0. (4.4) 

It turns out that the vector function {A2, B2) also satisfies a certain symmetry condition, as 
stated in the following proposition. 

Proposition 4.1 Let qx(0,y) be a sufficiently smooth function. Then the vector solution of 
the linear Volterra integral equation J^.iip satisfies the following symmetry conditions ( where 
we do not indicate the explicit dependence of A2, B2 on y): 

A2{\) = J— ^[^2(A) - i^(A)B2(A) + F(A)e-(^)(2'-^)i?2(-A) - F(A)2e-(^)(^-^)A2(-A)], 
S2(i) = ^_\^^xy [B2{\) - F{\)A2{\) + F(A)e"(^)(^-^)^2(-A) - F{\f&'^^^'^y-'^^ B2{-X% 



14 



M.,A)=( ), 0<.<L, AeC. (4.7) 



0<y<L, XeC, (4.5) 

where the function F{X) is given by 

, , 1 - d , ^ 

F{X)=i-^^tiin-. (4.6) 

Proof: Let the 2x2 matrix valued function $2(2/1 ^) be defined by 

A2iy,X) B2{y,-X) 
B2{y,X) A2{y,-X) 

Then $2 satisfies the ODE 

($2)y + ^[fT3, $2] = iQ{0, y, -A)$2, 0<y<L, 

<^2{L,X)^I, (4.8) 

where Q{x, y, X) is defined in (|2.5p . and (7(0, y) = d. 
Letting 

'i>2(y,A)-02(y,A)e^^-'(^-^) (4.9) 
it follows that (j)2 satisfies the ODE 

i4'2)v = VC^2. (4.10) 

(/)2(L,A)=/, 0<y<L, i?e(A)>0, 



where 



1/ -(A+^) -g.(0,y) 



i sin 



We seek a non singular matrix i?(A), independent of y, such that 



V{y,j) = R{X)V{y,X}R{Xr\ (4.12) 



It can be verified that such a matrix is given by 

1 

-FiX) 1 



m=[ I,, ^^^^ ), (4.13) 



where F is defined by (|4.6p . 

Replacing in equation (|4.10p A by i, and using (|4.12p . we find the following equation: 



Rixr'Mv, ^)) = V{y, A) (^R{x)-'My, ^)) , 



hence 

i?(A)-V2(2/,^) = <^2(y,A)C(A), 
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where C is a y- independent matrix. Using the second of equations (|4.8p . it follows that 
C = R~^, and therefore 

This equation and equation (|4.9p imply 

$2(y, j) = RWMy, A) (e-(^)^(^-^)i?(A)-i) . (4.14) 

The first column vector of this equation implies ()4.5p . 
QED 

Remark 4.1 Recalling that a2(A) = yl2(0, A), and 62(A) = -62(0, A), equations imme- 
diately imply the following important relations: 

a2{j) = Y^:^[«2(A) - F(A)fe2(A) + ^^(A)e--(^)H2(-A) - ^^(A)2e--(^)^a2(-A)], 

b2{\ = Y-r^M>^) F{\)a2{\) + ^^(A)e--(^)^a2(-A) - Fi\fe-'^WH,i-X)], 
Im{X) > 0. (4.15) 
In summary, the basic equations characterizing the spectral functions are: 

(a) the symmetry relations (|4.4[) and (|4.15p : 

(b) the global relations (^T^ and (I2221); 

(c) the conditions of unit determinant. 

It turns out that, using these equations, it is possible to provide an explicit characterization 
of all the spectral functions in terms of the given constant d. 

Proposition 4.2 Assume that the functions {aj{X), bj(X)}, j = 1, 2, 3 satisfy the symmetry 
relations ^4-4^ o-i^d, ^.15^ , the global relations h2.21]) and \2.22\) and the "unit determinant" 
conditions 



a,(A)a,(-A)-6,(A)5,(-A) = l, 


= 1,2,3. 


(4.16) 


Then the following relations are valid: 






ai(A)6i(-A)-ai(-A)6i(A) = G(A), 


A G R, 


(4.17) 


a3(A)63(-A) - a3(-A)63(A) = -G(A), 


A e M, 


(4.18) 


a2(A) = ai(-A)a3(A) - c— (^)^6i(A)63(-A), 


A e C, 


(4.19) 


62(A) = ai(-A)53(A) - e— (^)^6i(A)a3(-A), 


A e C. 


(4.20) 


where 

i(l-A2)e"(^)^+e-"(^)^-2 


d 

tan-. 


(4.21) 
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Remark 4.2 The two relations J^.Jffl j and ^.20^ are a direct consequence of the global 
relation and of the conditions of unit determinant. On the other hand, equations {jj-l'T^ and 
depend on the particular symmetry properties. 

Proof: for simplicity, we will use the notations 

/ = /(A), /-/(-A). (4.22) 

Replacing A with —A in (|2.2ip and solving the resulting equation and equation (|2.22p for 
a2(A) and &2(A) we find equations (|4?T9)l and ((420)1 . 

It can be verified directly that if {02,62} are defined by equations (|4.19p and (|4.20p and 
{aj,6j}, j = 1,3 satisfy the determinant condition (|4.16p . then {02,62} also satisfies the 
determinant condition. 

Replacing in the global relations (|2.2ip and (|2.22p A by — and using in the resulting 
equations the symmetry relations (j4.4p . we find 



ai(A) = a2(^)a3(A) - h2{\)h{\). (4.23) 

6i(A)e-(^)^ = a2(-i)63(A) - b2{-\)a^{X) . (4.24) 

Replacing in these equations 02 (±j) and 62 (ij) by the right hand side of the symmetry 
relation (|4.15p . as well as by the right hand side of the equation obtained from equations 
(|4.15p under the transformation A —A, after extensive simplifications we find the following 
equations, valid for A G R: 

{l-F^)ai 2 ^2, ai-F6i , ^ . . 

b. Fa, - - + ®' ^'-'^^ 

il-F^)b, ,2, bi- Fa, , ^ 

a, - Fb, = - ^^)^7^ + ®' ^'■2^) 

where 

e = c-"(^)^(63a3 - 6303) + e-'^^^^^F, A G R. 

Hence 

(1 - F'^)ai + A(ai - Fb,) (1 - F'^)bi + A(6i - Fa, 



bl — Fai ai — Fbi 



e, A 63 - 4. (4.27) 



Multiplying the numerator and denominator of the first fraction in (I4.27P by F and adding 
the resulting expression to the second fraction, as well as multiplying the numerator and 
denominator of the second fraction in (|4.27p by F and adding the resulting expression to 
the first fraction, we find 

Fai+61+ASi ai+Fbi + Aai 

- - fc). (4.2«j 

ai bl 

Using equation (|4.16p with j — I, the left hand side of (I4.28P implies 

7,2 2 7,2 2 
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Then, the second equation of (I4.27P imphes 

oi ^ a,ibl -al) ^ ^-.(A)L(^^^^ _ ^ ^-^WLp _ p_ (4,29) 

Oi Oi 

The left hand side of (|4.29p . using (|4.16p with j ~ I, simphfies as follows: 

^ + aibi - ^(1 + bibi) = aibi - aibi. 
bi bi 

Thus, equation l|4.29p becomes 

fciai - hai = e-'^(^)^(63a3 - 6303) + e-'^'-^^^F - F. (4.30) 

Replacing in this equation A by —A yields 

- {biai - hai) = -e"(^)^(63a3 - 6303) + e'^^^^^F - F. (4.31) 

Equations (|4.30p and (|4.3ip . taking into account the definition (|4.6p of F, yield equations 

(jilT)) and Km . 

QED 

Remark 4.3 The determinant condition (I4.16P with j ~ 1 and equation (|4.17p imply 

[(ai(A)2-6i(A)2)] [(ai(-A)2-6i(-A)2)] =1-G(A)2, A e M. (4.32) 
Indeed, equation (j4.16p with j = 1 implies the identity 

[{af - bl)] \{al ^ bj)] = I - {aih - aibif , A G M. (4.33) 

Then equation l|iT7|) implies 

Equation (|4.32p defines the jump relation of a scalar RH problem for the sectional analytic 
functions defined by 

{(ai(A)2 - 6i(A)2, A e C+ ai{-Xf - bi{-\)\ A e C-}. 

Taking into consideration that ai(A) 7^ bi{X) for A £ C+ (otherwise equation (|4.16p with 
= 1 is violated) it follows that the above Riemann-Hilbert problem has a unique solution 

ai(A)^ -6i(A)2 = /i(A), AeC+, (4.34) 

where 



MA)=e^(^), H{X)^— \n[l-G'{\')]-^, A e C \ K. (4.35) 



1 / ln[l-G^(A')]-^., 
2711 A' - a 

Using the fact that G(A) is an odd function, it follows that H{X) is also an odd function, 
hence h{—X) = e~^^'^\ This implies that the function h{X) defined by (I4.35P satisfies the 
jump condition (|4.32p . 
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Remark 4.4 Equation (|4.16p with j = 1, and equations (|4.17p and (|4.34p imply 

«i(-A)-^(ai(A)+G(A)6i(A)), 6i(-A) = ^ (61(A) + G(A)ai(A)) , A e M. 

(4.36) 

Indeed, equation (|4.17p yields 

61 = — (G + ai6i). 

ai 

Replacing bi in equation (|4.16p with j = 1 by the above expression, and making use of 
(|4.34p . we find the first of equations (|4.36p . The second of equations (|4.36p can be obtained 
in a similar way by eliminating oi instead of bi. 

Remark 4.5 The equations satisfied by 03 and 63 can be obtained from equations (|4.36p 
by replacing G(A) by G(— A). Hence 

a3(-A) = ^(a3(A)-G(A)63(A)), 63(-A) = ^ (63(A) - G(A)a3(A)) , A e M, 

(4.37) 

where G is given by (|4.2ip and h{X) is given by (|4.35p . 

Remark 4.6 The function G is an entire function, thus each of equations (|4.36p defines 
the jump condition of a scalar RH problem. However, it will be shown in section 5 that 
equations (|4.36p and (|4.37p are sufficient to determine the jump matrix (|2.25p . 

Remark 4.7 Equations (|4.17p - (|4.20p imply the following identity: 

e"(^)^[ai(A)2~6i(A)2]+e-'^(^)^[ai(-A)2-6i(-A)2] = (c"(^)^+e-'^(^)■^)(l-F2)+2i^^ XeR. 

(4.38) 

Indeed, equations imply 

gC.(A)L(^2 _ ^2) ^ g^(A)L^2(^2 _ ^2) _ g-^(A)L^2(^2 _ ^2) _ 2^161(0363 - 6363) (4.39) 

Replacing in this equation A by —A, adding the resulting equation to equation (|4.39p and 
using equation (|4.30p we find 

e^^^^'^{al-bl)+e-'^'^^^'^{al-bl) = {e'^<^^^'^+e'^^^^^^ 

(4.40) 

Using (|4.32p . the right hand side of (|4.40p equals the following expression: 

(g^(A)L ^ g-c.(A)L) _ (^^^^ _ J(e-W^ + e-"(^)^)(ai6i - a^bi) - 2(0363 - 03^3) • 

Using equations (|4.17p and (|4.18p the last expression becomes the right hand side of (|4.38p . 



5 Spectral theory in the Hnearisable case 

In the case of the linearisable boundary conditions (|1.2p . it is possible to express q{x, y) in 
terms of the solution of a RH problem whose jump matrices are computed explicitly in terms 
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of the given constant d. Indeed, recall that the jump matrices of the basic RH problem of 
section [2^ are defined as follows: 



J 



tt/2 _ 



1 /(A)e-''(^'^'^) 



j3tt/2 _ 




/(A) 



fe2(-A) 
ai(A)a3(A) ' 



R{X) 



ai(-A) ^ 



1 



R{X) = 



a2(A) 



ai(-A)a3(A)' 



and 



Equation (|4T9| and (|420| imply that 



i?(A) 1 -e" 



ai 



03 ai 03 

^and^ 

at A and at —A. Equations (|4.35p and (|4.34p imply that these rations are given by 



(5.1) 



(5.2) 



Thus in the linearisable case, the jump matrices involve only the rations and evaluated 



as 



h a^h ' 



h_G bi_ 
fli h aih 

G 



(5.3) 



Hence the jump matrices depend on the known function ^ as well as on the unknown 

Using the fact that these unknown functions are bounded and 



functions 



ai h 



and 



analytic in , it is possible to formulate a RH problem, equivalent to the basic one defined 
by (j5.ip , in terms of the known function ^ only. This new RH problem is therefore defined 
by the following jump matrices: 



1 /(A)e-^(^'2'^^) 



R{\) 



j37r/2 



^G_f.e{x.,v,\) 




7(A) = -^(1+e- 



.(A) 



h 



hh 



o(\)L 



and 



(5.4) 



Theorem 5.1 Let q{x,y) satisfy equation il.l}) and the boundary conditions lil.2\) . 
Then q{x, y) is given by equations i2.36\ )- (2.S7^ with ^! replaced by ^/ , where 4* is the solution 
of the Riemann-Hilbert problem i2. 23]) with the jump matrix J replaced by the matrix J 
defined as follows: 



J{x,yA) = J°'ix,y,X), ifarg{X) = a, a = 0, ^, tt, 



(5.5) 
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wher 



1 _ G'(A) -u^(\)L _G(A) e(x,a,A) 
h(X)h(-X)^ TtJ)^ 



-uj(\)L G(A) e(x.v.\) 



/l „G(A2n+e^(A)L-)g-e(.,,,A) N / 1 

V 1 / ' V ^(l + e~-W^)e«(-.«.^) 1 

and 

J- = J3V2(J0)-1 jV2^ (5 6) 

where 9{x, y, A) is given by i2.26]} . while G{X) and h{X) are defined in terms of the given 
constant d by equations {4-21^ and \4-.35^ . 

This Riemann-Hilbert problem is regular and, ij d ^ R, it has a unique solution. 

Proof: The solution of the new RH problem satisfies the jump relation 

y, A) = y, A) J(x, A), A e M U iR, (5.7) 

where the jump matrix J is given by equations (j5.5p and (j5.6p (the latter are identical to 
equations (|5.4p ). Let ^'j and j = 1, ...,4, denote and 5* in the j-th quadrant of the 
complex A plane, respectively. We seek matrices Aj, j — 1,...,4 which are analytic and 
bounded in the j-th quadrant of the complex A plane, respectively, and such that 



^j{x,y,\) = ^j{x,y,\)A,{x,y,\), {j - l)^ < arg{\) < j^, j = l,...,4. (5.8) 



Let 



1 ai(A)e 



Ai= \ I , < argiX) < J (5.9) 



1 

1 a2(A)e-^(^'2''^) 



2 



TT 



A3 = I I , < "^5(A) < ^ (5.11) 

c,3(A)e'^(-.a,A) 1 J ^ 

A4= \ 1 ' ^ < «''5(A) < 27r. (5.12) 

a^{\y(^-v^) 1 y ^ 

Equations (|2.23|) and (|5.7|) imply the following relations: 

J'^/^Aa = AiJ""/^, J^^'/^Ai ^ AsJ^""/^, J°A4 = AJ°. (5.13) 
The first two equations (|5.13p imply 

a2(A) +/(A) = ai(A) +/(A), a4(A) + /(-A) a3(A) + /(-A). (5.14) 
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Using the definitions of J(A) and J(A) (see equations (|5.2p and (|5.4p ') as well as (|5.3p . equa- 
tions (|5.14p become 

a2 H r — e ^ ' — - = ai, Q!4 H ^ — e ^ ^ — - — a^. 

The simplest solution of these equations that satisfy the requirement that the functions 
ocj{X) are bounded and analytic in the j-th quadrant of the complex A-plane, j = 1, ...,4, 
are the following: 

a^n aih 

as = — a4 = — ^e ^ ' . (5.15) 
0,3/1 aih 

The third of equations (|5.13p yields the following relations: 

as ft ai h 

R + a4— ^ R + aie-^^^^^^. (5.16) 
a3 /i 

Using the definitions of i?(A), i?(A) (equations (|5.ip . (|5.4p ) the definitions of , j = 1, ..,4 
(equations (|5.15p ) and and equations (|5.3p it can be verified that equations (|5.16p are satisfied 
identically. 

The eigenfunctions ^'j, j = 1, ...,4 are given explicitly in terms of the eigenfunctions by 
equations (|5.8p . This yields the following expressions: 

03/1 03 aih ai 

(5.17) 

§3 ^ (^vff + iLvl/f )e«(-.^.^), vl/f )), ^4 ^ (Ivl/W + _^e-(^)^vl/f )ge(.,,,A)^ ^(34)^^ 
03 0,3/1 Ol Oi/l 

The above equations can be simplified as follows: 

(5.18) 

§3 = (^vi/f) + ^vi,We«(-.s^.>), vff )), §4 ^ (^vi/(4) + bie-(^)^ ^(2)^g(,,^,^)^^(34)^_ 
h h h h 

Indeed, using equation (|2.16p we find 

_fe^^(i2)^_e(x,,,A) ^ 1^(1) ^ _^(a3M/(3) +53*«e'^(-'^'^))e-''(-'^'^) + 

03/1 03 03/1 J J ^3 
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Using a§ — 6| = /i the above expression is equal to the regular function 

A similar computation yields the result for the other eigenfunctions. 
Thus the above Riemann-Hilbert problem ()5.7p is regular. 

We now prove that the Riemann-Hilbert problem defined by (|5.7[) is uniquely solvable. It 
can be verified that when d G M, then h{\) = h(X). In this case, the jump matrices J'"-' 
satisfy the following conditions: the matrices are Schwarz invariant on the imaginary axis 
and have zero real part on the real axis of the complex A plane. Under these assumptions, 
it follows from general results ( see e.g. [TUl that the so-called "vanishing lemma" 

holds. This guarantees the existence of a unique solution. 
QED 

6 Conclusions 

We have studied boundary value problems for the elliptic sine-Gordon posed on a semistrip. 
In particular we have shown that if the prescribed boundary conditions are zero along the 
unbounded sides of the semistrip and constant on the bounded side, then it is possible to 
obtain the solution in terms of a Riemann-Hilbert problem which is uniquely defined in 
terms of the width L of the semistrip and the constant value d of the solution along the 
a; = boundary. Indeed, the "jump matrices" of this Riemann-Hilbert problem are defined 
in terms of the two functions G(A) and h{X) which are explicitly defined, by equations (|4.2ip 
and (j4.35p respectively, in terms of L and d. Furthermore, this Riemann-Hilbert problem 
has a unique solution, see theorem 5.1. 
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